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We study the effect of left-right symmetry on the distribution of the parametric derivative of the 
dimensionless conductance T with respect to an external parameter X, dT/dX, of ballistic chaotic 
cavities with two leads, each supporting N propagating modes. We show that T and dT/dX are 
linearly uncorrelated for any N. For N = 1 we calculate the distribution of dT/dX in the presence 
and absence of time-reversal invariance. In both cases, it has a logarithmic singularity at zero 
derivative and algebraic tails with an exponent different from the one of the asymmetric case. We 
also obtain explicit analytical results for the mean and variance of the distribution of dT/dX for 
arbitrary N. Numerical simulations are performed for N — 5 and 10 to show that the distribution 
P(dT/dX) tends towards a Gaussian one when N increases. 

PACS numbers: 05.45.-a, 42.25.Bs, 73.23.-b 



I. INTRODUCTION 

In recent years there has been great interest in coher- 
ent quantum transport in mesoscopic devices in both ex- 
periment and theory (for a review see Refi). In quan- 
tum systems whose classical dynamics is chaotic, phase 
coherence and quantum interference give rise to samplc- 
to-sample fluctuations in most transport properties with 
respect to small perturbations in incident energy, cavity 
shape, or an external parameter such as an applied mag- 
netic field. Therefore an statistical analysis by random- 
matrix theory (RMT) becomes applicable; the details of 
the problem being irrelevant emphasis must be put on 
universal aspects such as the symmetries of the system. 
Then, time-reversal invariance (TRI) and spin-rotational 
symmetries^, or the spatial symmetries present in the 
system characterize their universal statistics 3 ' 4 ' 5,6 . 

Some experiments concerned with the parametric de- 
pendence of the conductance study ballistic quantum 
dots connected by leads with few propagating modes 
to electron reservoirs, as happens in semiconductor 
systems2*2i2iiSiii. The conductance is measured as a 
function of the magnetic field and the shape of the quan- 
tum dot (controlled by gate voltages). The derivative of 
the conductance with respect to an external parameter, 
which is analogous to the level velocitjiiSiiiiiii^ of the 
dot, is a very important quantity that gives us the re- 
sponse of a mesoscopic sample to a small perturbation. 
Only the asymmetrical situation has been considered in 
the literature using an RMT approachiiiSiiiiS* 1 ^. 

Other methods to verify RMT predictions are given by 
wave scattering experiments, such as microwaves22i2ii22 
and acoustic resonators^ 3 -. In the microwave technique, 
for example, the conductance is measured as a function 
of the frequency of incidence, applied magnetic field, and 
cavity shape. 

In these wave scattering experiments the external pa- 
rameters are easier to control. The symmetric case can 
be constructed to analyze the strong effects of spatial 



symmetries on the conductance, predicted by RMT^. 
Experiments have been conducted on left-right symmet- 
ric cavities^ 2 . However, no theoretical results exist for the 
distribution of the conductance velocities. The purpose 
of this paper is to study the effect of the spatial left-right 
symmetry in the distribution of the parametric derivative 
of the conductance. Thus, we consider a chaotic left- 
right symmetric cavity connected to two symmetrically 
located leads for two different symmetries, with and with- 
out TRI, in the absence of an applied magnetic field 24 . 
Other spatial symmetries^ are not considered. 

When TRI is broken by a magnetic field the problem 
of left-right symmetric cavities is mapped to the one of 
asymmetric cavities in the presence of TRI. That problem 
was studied by Brouwer et a& using an RMT method to 
calculate analytically the joint distribution of the dimen- 
sionless conductance T and its velocities dT/ dV, dT/dX, 
with respect to the gate voltage V and a external param- 
eter X (typically a magnetic field), for a quantum dot 
with two single- mode point contacts. The statistical dis- 
tributions P(dT/dV), P{dT/dX) show algebraic tails, 
and in the absence (presence) of TRI they have a cusp 
(divergence) at zero velocity, and the second moment is 
finite (infinite). RMT predictions have been verified by 
experimental results^. 

In the present work we perform analytical calculations 
for the one channel case to reduce to quadratures the 
distributions of dT/dE, and the derivative dT/dX with 
respect to an external parameter X (shape deformation 
parameter in microwave cavities), using a scattering ma- 
trix formalism. By direct numerical integration we obtain 
P (dT/dE), P(dT/dX). Analytically, we obtain their 
mean and variance and show also that T, dT/ dE, dT/ dX 
are linearly uncorrelated for any number of channels. Be- 
sides, we obtain P(dT/dE), P(dT/dX) by numerical 
simulations using an scattering matrix. 

The formal elements such as the scattering matrix for 
the symmetric case, its energy and parametric deriva- 
tives are introduced in the next section, as well as the 
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FIG. 1: A ballistic chaotic left- right symmetric cavity con- 
nected to two leads supporting one channel each one. 



joint distribution of them, for an arbitrary number of 
propagating modes TV in the leads. The TV — 1 case is 
the subject of Sect. 1111 Al for both in the presence and 
absence of TRI. Sect. IIII Bl is dedicated to an arbitrary 
TV case. Finally, the conclusions are presented in Sect. 



CHAOTIC SCATTERING BY REFLECTION 
SYMMETRIC CAVITIES 



II. 



The scattering problem of a ballistic cavity with left- 
right (LR) symmetry connected to two leads, each with 
TV transverse propagating modes is described by the 
2TV x 2TV scattering matrix S, which in the stationary 
case relates the amplitudes of the outgoing to the incom- 
ing plane waves 2 ^. For a system with reflection symmetry 
(see Fig. ^| the S matrix is block diagonal in a basis of 
definite parity with respect to reflections. The general 
structure for 



S = 



r t 
t r 



(1) 



where r, t are the TV x TV reflection and transmission 
matrices. 

The transmission coefficient, or spinless dimensionless 
conductance, is obtained from the S matrix 

r = tr(ttt) ; (2) 

it is proportional to the conductance of the cavity, 

G = (e 2 /h)T. (3) 

Matrices with the structure can be brought to the 
block-diagonal form 



S ~ R ° ( o 1 S 2 ) R °> 

where Rq is the rotation matrix 

p _ _L ( I* ijv\ 

l n denotes the n x n unit matrix and 

r = (S 1 + S 2 )/2 
t = (Si - S 2 )/2. 



(4) 



(5) 



(6) 
(7) 



Here, Si = r + 1, S 2 = r — t are the most general TV x TV 
scattering matrices. 

In the absence of any symmetry the condition of flux 
conservation implies unitarity for Sj (j = 1,2): SjSj — 
In- In the Dyson's scheme^ this case is called "unitary" 
and it is designated by (3 = 2. In addition, in the presence 
of time reversal invariance Sj is symmetric: Sj = Sj '. 
This is the "orthogonal" case, designated as = 1. It 
means that two independent N x N scattering matrices 
Si, S 2 unitary (and symmetric) for (3 = 2(1) generate 
the most general 27V x 27V unitary (and symmetric) S 
matrix with the structure (Q: 



SS^ = 1 2N , 
S — s • 



for (3 
for P 
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(8) 
(9) 



When the classical dynamics of the system is chaotic, 
most transport properties are sample specific and a sta- 
tistical analysis of the quantum-mechanical problem is 
called for. That analysis is performed by the construction 
of ensembles of physical systems, described mathemati- 
cally by ensembles of S matrices distributed according 
to a probability law. The starting point is a uniform 
distribution where the scattering matrix is a member of 
one of the circular ensembles, defined by the invariant 
measure which is the precise formulation of the intuitive 
notion of equal a priori probabilities in the space of scat- 
tering matrices. In our case Si and 5*2 are independent 
of each other, statistically uncorrelated, and distributed 
according to the invariant measure defined by 



dfi^(Sj) = dfx^(U SjV ), (j = 1,2). 



(10) 



Here, Uq, Vq are arbitrary but fixed unitary matrices for 
(3 = 2, while V = Uq for ft = 1. Equation JJ3 defines 
the circular unitary (orthogonal) ensemble, CUE (COE), 
for 13 = 2 {(3 = 1)21. 

Then, the invariant measure for S with the structure 
||TJ can be written as^ 



dp^\S) = d^{Si)d^\S 2 ). 



(11) 



S matrices of the form given by Eq. QJ, which sat- 
isfy Eq. (JSJ) are appropriate for systems with reflection 
symmetry in the absence of TRI. With the additional 
condition |J3J| it is appropriate for LR-systems in the pres- 
ence of TRI. However, when TRI is broken by a uniform 
magnetic field, the problem of LR-symmctric cavities is 
mapped to one of asymmetric cavities* with (3 = 1 with 
t replaced by r. We are interested in the case when 
TRI is not broken by a magnetic fieldSi. Other spatial 
symmetries^ will not be considered here. 

The parametric derivative of Sj (J = 1, 2) with respect 
to the energy of incidence can be defined in terms of a 
symmetrized form of the Wigner-Smith time delay ma- 
trix, whose eigenvalues are also the proper delay times. 
The derivative of S with respect to another external pa- 
rameter X is defined in similar way. We have 28 



dE 



i sj Qj E sj , 



dX 



1 S) /2 Q 



x Sj 



1/2 
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Qj E , Qj x are N x N Hermitian matrices for P = 2, real 
symmetric for [3—1. The eigenvalues of Qj E are 
times the proper delay times. 

For classically chaotic cavities the joint distribution of 
Sj > Qj e i Qj v i^"^ 



Pf ] {Sj,Q jE ,Q jx )cx(detQj E 



x exp < —p tr 



'3E 



-2/3Af-3(l-/3/2) 
2" 



2A 



Q?\e Q 



JX 



, (13) 



where A is the mean level spacing, and Xp is a typi- 
cal scale for X which is not a universal quantity. Sj is 
independent of Qj E and Qj x and it is uniformly dis- 
tributed in the space of scattering matrices, i.e. ac- 
cording to the invariant measure dfJt^'(Sj). Qj x 1S dis- 
tributed as a Gaussian with a width set by Qj E - The 
rates {xj — 1/rj } (n = 1, ... ,N), where the dimen- 
sionless time delay tj is A/2tt times the nth eigenvalue 
of Qj E , are distributed according to the Laguerre ensem- 
ble, namely^ 



Pf> far}) « n 



J a 



°3b\ 



n 



(14) 



a<b 



For N = 1 it reproduces the result of Refsi 16 i 18 i 29 . 

Let denote by q any of the two parameters E or X. 
A convenient parametrization for Sj and dSj/dq (q — 
E, X) is 



Si 



dS 



3 =iU j Q jq V j , 



(15) 



where Uj , Vj are the most general Nx N unitary matrices 
for = 2, while Vj = Uj for p = 1. 
We parametrize Qj X and as2£ 



0. 



J.Y 



27T if i 
A J J 



(16) 
(17) 



where Kj is & N x N Hermitian matrix for p = 2, real 
symmetric for P = 1, that have a Gaussian distribution 
with zero mean and variance 

((i ^" )ab( ^ u = { mtt + s ac s bd ) pZl > ( 18 ) 

as can be seen by substitution of Eqs. I|16|) . (|17|) into Eq. 
()13|) . 'J is a A~ x N matrix, complex in the unitary case, 

real in the orthogonal one. We define Qj = tyj vj/7 1 
and diagonalize it, 



Qj = WjTjW] 



(19) 



The reciprocals of the eigenvalues {Tj } of fj are dis- 
tributed according to Eq. I|14|). The matrix of eigenvec- 
tors, W, is uniformly distributed in the unitary (orthog- 
onal) group for P — 2 (/3 = 1). 



III. PARAMETRIC CONDUCTANCE 
VELOCITY DISTRIBUTIONS 

A. The iV = 1 case 

In this section we calculate the distributions of the par- 
tial derivatives of the dimensionless conductance, namely 
dT/dE and dT/dX, when the leads support one propa- 
gating mode. 



For N 



1 the S matrix is 2 x 2 and Si 



(j = 1,2) is uniformly distributed in the unit circle, i.e. 
dn^\Sj) = d6j /2tt; S is distributed as Eq. Q) 



d£W(S) 

The transmission coefficient is 
1 



(201 d(9 2 



Hi 



whose probability distribution is given by2> 
wp(T) 



n^T(\-T) 



(20) 



(21) 



(22) 



for any p. 

The partial derivative with respect to the external pa- 
rameter q (q = E, X) is 



dT 



1 



d6 x 802 



from (|21|l we rewrite the last equation as 
1 dT ( 50, 96» 2 



cr^T(l - T) #9 



(23) 



= - > ( 24 ) 



<90_i 
9g 



where <r = ±1 depending on the phases of the 1 X 1 Si, S2 
matrices, 9\, 62- A similar expression to Eq. I|24|) holds 
for asymmetric cavities^ where dT/dE depends on T 
but the ratio (dT/dE)/y/T(l - T) is independent of T. 
However, T is linearly uncorrelated with dT/dq (q = E, 
X) (see Sect. EH). 
We write Eq. (J23J) as 



dT 



= — sin 

2 V 



2). 



where we have defined 



<96>i «96» 2 



(25) 



(26) 



From Eq. (|12J) . dOj/dE = Tj . In similar way we denote 
by = dOj/dX. Then, 



ze = n - r 2 , 2x = Pi - P2- 



(27) 



The sets (6*i, Ti, pi), (0 2 , T2, P2) are statistically indepen- 
dent, while the joint distribution of 9j, Tj, pj (J = 1,2) 
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is [Eq. 



P 



,pj)= C P T 



-3-/3/2 



exp 



ap a p x p P) 
r, + 4/3 rf 



_ (28) 

where ap = 0n/A, Cp = a (2+(3/2) Xp/2y/0lrT(l + 0/2), 
and T(x) is the gamma function. 

The distribution of dT/dq (q = E, X) is 



2tt 



2tt 



x J 



2tt 



dz q P'p{z q ) 



which can be written as, 

Pp(dT/dq) = - 
or, equivalently, 
Pp(dT/dq) = 



P' 



2 <9T\ d0 
sin 6* <9g y | sin6>| ' 



(29) 



(30) 



n J2\dT/dq\ - 4(<9T/«9g) 2 



dz q , (31) 



where 



2w d9i 
2tt 

2 " 



'0 

x S[z, 



2tt 

, - (6 



dr 2 



(/3), 



(32) 



with £ = r or p. Using Eq. H28[l we see that the inte- 
grations with respect to 6\, 9 2 are equal to 1. For the 
particular case of Pp(dT/dE) the p variables are easily 
integrated, as well as one of the r variables, say T2, due 
to the delta function factor, thus only remaining to be 
done the integral with respect to t%. After this is done 
and using Eq. l(37)|) we get that 



Pp(dT/dE) = A t 



dT 


-(2+/3) , 


dE 


Jo 



dv- 



lll^+//3M>f§)-sin 



d9 (sm6)( 1+ W 

2+/3 



/MM, f!) \m\ £ + fp( v > e >m) 



1+/3/2 • 



(33) 



where fp(v,9,dT/dE) = {[{dT / dE)v / ap] 2 + sin 2 0} 1 / 2 
and ^ = aJj 1+/3) /8/?7r[r(l + /5/2)] 2 . 

For P(dT/dX) one of the p's, say p%, is easily inte- 
grated due to the delta function factor, giving a Gaussian 
integral for p\ which is also easily obtained. The result 
is substituted in Eq. I|31|l , the order of integration is in- 
terchanged to first perform the Zx integration, using an 
integral representation of the zero order modified Bessel 
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FIG. 2: Distribution of (a) AdT/dE and (b) XpdT/dX for 
the symmetries = 1,2, obtained by numerical integrations 
of Eqs. |SSI and 



function Kq(x). The remaining integrals, with respect to 
Ti , T2 , are transformed to new variables to obtain that 



Pp(dT/dX) = Bp f dx 
Jo 



(1-x) 



(1+/3/2) 



dy- 



x (2+p/2) 
(1 - y)(l+/>/3) 



-a l 3(l~x)/x 



y(2+/3/2) ^ _ x y y 2 + x 2(i _ y y 

2 



(!— y)/y 



x exp 



' dT 

-9p{x,y) ( ^ 



Ko 



dTY 
9ff(x,y) [ — J 



(34) 



where Bp = 2Xpaf +P) /ny/0^[T{l + 0/2)] 2 and 



9p{x,y) 



a 2 X 2 



{l-x) 2 (l-yf 



20 (l-x)V + x 2 (l~y) 2 



(35) 



In Fig. |21 we show the result of the numerical integra- 
tion of Eqs. H33fl an d 134|l . We observe a slight depen- 
dence on 0, while the T-distribution is independent of it 
[see Eq. (|2*2*)) ]. In our case, a symmetric one, the depen- 
dence on is not as clearly marked as that of the asym- 
metric case. Also, in contrast to the asymmetric casein, 
Pp(dT / dX) diverges logarithmically at zero transmission 
velocity not only for — 1 but also for = 2. 

For large values of \dT/dq\ the tails are algebraic, 
namely 



Pp{dT/dq) <x \8T/dq\ 



-2-/3/2 



(36) 
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which again is in contrast with the asymmetric case 
where the exponent of the algrabraic tails^S is —2 — (3. 



B. Arbitrary N 

For arbitrary N we are not able to calculate analyti- 
cally the whole distribution of dT/dq (q = E, X). In- 
stead, we calculate the first two moments of the distribu- 
tion analytically for any N and complement our results 
with a numerical simulation of Pp(dT/dq) for N = 5, 
10. Also, we calculate the covariance between the energy 
and parametric derivatives, as well as the correlation of 
each one with the transmission coefficient. The explicit 
analytical calculations are in the Appendix. 

Due to Eq. (0J the transmission coefficient depends 
on the interference term SiS^ (see Eq. lAlfl . i.e. linearly 
on S± and S2, that are independent from each other, 
but are equally distributed random matrices. As a re- 
sult we have (see the Appendix) (T{dT / dq)) = and 
((dT/dE)(dT/dX)) = for all N. This means that 
T, dT/dE, dT/dX are linearly uncorrelated variables. 
However, in general they are correlated (see the Ap- 
pendix [SJ as happens in the asymmetric casein. 

We find (Appendix) that Pp(dT/dq) has a zero mean, 
i.e. (dT/dq) = 0, as expected because P/3(dT/dq) is a 
even function of dT/dq (see Fig. 0). For N = 1, which is 
the case of Sect. IIII Al ((dT/dq) 2 ) diverges in both cases 
((3 = 1,2). For N > 1 the variance of dT/dX is given by 



N=10 



(dT/dX)' 



for (3 = 1, and 



(dT/dxf 



for (3 = 2; for both cases 



1 2(N 2 + N + 2) 
Xj (N- 2)(N+ l) 2 



1 N 2 + 1 
Y 2 N(N 2 - 1) 



(dT/dE)' 



A 2 N 



(dT/dX)' 



(37) 



(38) 



(39) 



We see that ((dT/dq) 2 ) diverges only for j3 = 1 when 
N = 2. 

For large N cases we resort to numerical simulations. 
The results for N = 5, 10 are presented in Fig. [21 for 
dT/dX, and in Fig. 0]for dT/dE where we also include 
the Gaussian distributions with variances obtained from 
our results given by Eqs. (33 and jSHJ for Pp(dT/dX), 
and the corresponding ones for Pp(dT / 'dE). We observe 
that the distribution tends towards a Gaussian as N is 
increased as expected (see Ref^ and references there in). 

As mentioned in a previous publication—, for chaotic 
cavities we implement numerical simulations by generat- 
ing 10 4 scattering matrices for each Sj (j = 1, 2) given by 
the Heidelberg approach^ with M = 100 poles. Details 
can be found in Refi^i. 
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FIG. 3: Distribution of dT/dX for N = 5, 10. When N 
increases the distribution tends to a Gaussian distribution 
with variance given by Eq. 11371 for (3 = 1 and Eq. 1381 for 

13 = 2. 
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FIG. 4: As in Fig. El the distribution of dT/dE tends to a 
Gausssian distribution with variance given by Eq. 1391 . 



IV. CONCLUSIONS 

We have studied the effect of spatial reflection sym- 
metry on the distribution of the parametric derivative, 
dT/dX, of the dimensionless conductance T of a chaotic 
cavity, with respect to an external parameter X. In par- 
ticular we considered a left-right ballistic chaotic that 
has two leads symmetrically placed, each one supporting 
iV propagating modes, in the presence (absence) (3 — 1 
((3 = 2) of time-reversal invariance (TRI). The external 
parameter could be, for instance, the shape of the cav- 
ity in the microwave technique, where our results can be 
tested^. 

For the particular case N — 1, we calculated the distri- 
bution Pp(dT / dX) by reducing the problem to quadra- 
tures and performing a numerical integration. While the 
distribution of T does not depend on /3, the distribution 
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of dT/dX show a slight dependence on (3 that is not as 
clearly marked as in the asymmetric case. Also, in our 
case, Pp(dT / dX) diverges logarithmically at zero trans- 
mission velocity for any (3. This is in contrast with the 
asymmetric case where for (3 = 1 there is a logarithmic 
divergence and a cusp (3 = 2. We found that P/3(dT/dX) 
has algebraic tails with an exponent —2 — (3/2 different 
from that in the asymmetric case. 

We found that dT/dq (q = E, X) depends on T but 
the ratio (dT/dq)/ y/T(l - T) is T-independent as hap- 
pens in the asymmetric case. We showed that T and 
dT/dq are linearly uncorrelated for any N, although in 
general they are correlated. 

We performed numerical simulations for N = 5, 10 to 
obtain Pp(dT / dX) and compare with a Gaussian distri- 
butions with variances calculated with our analytical re- 
sult; we see that the distribution of the parametric deriva- 
tive of T tends towards a Gaussian distribution when N 
becomes large. 

Finally, we emphasize that our work is concerned with 
fully chaotic ballistic cavities. Our results are consistent 
with both experimental findings and simulations 22 . It 
would be also very interesting to single out the contri- 
bution of regular orbits (see^^S*^ and references there 
in) something that is beyond the scope of the present 
work, since it would require to modify the method here 
employed. 
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APPENDIX A: MEAN AND VARIANCE OF 

dT/dq, (q = E, X) 



From Eqs. @ and J7J| we see that 



N 



T = ^- \ E Re(SiU(s 2 *U 



(Al) 



.,6=1 



where we have used the unitarity condition for S\ and 
S2. Because Si and S2 are independent but equally dis- 
tributed random matrices with random phases, the aver- 
age of T is (T) = N/2 for both symmetries ((3 = 1 and 
2), reflecting the spatial symmetry of the cavity. 

Let q be one of the two parameter E or X in dimension- 
less units. For the same reason as in the last paragraph, 
the average of the g-derivative of T, 



dT 1 f p 



a, 6=1 



(Sl) a , 



d(si) ab , d(Si) ab 



dq 



+ 



-(S^al 



(A2) 



is (dT/dq) = 0. 

In similar way, after averaging the product of Eqs. 
(fAT) and ljA"2|l we get 



T 



,<9T 
Oq 



N 



Ke((SiU(SZ) a >v 



(S 



2/ab- 



d( S 2 )a>. 

dq 



(A3) 



The first term on the right hand side is real positive. For 
q = E, the second term of Eq. I|A3[) is pure imaginary as 
can be seen by substitution of Eqs. (|15|) . (|17l) and H19|) 
and averaging over the unitary matrices U2, V2 and W2 
for (3 = 2, or over U2 and W2 for (3=1. For q = X 
the parametrization Ijl5(l and Eq. Ijl6|l gives a result 
linear on the matrix K2 which has zero mean. There- 
fore, (T (dT / dq)) = 0. Following the same procedure we 
expect that (T m (dT/dq)) = with m an integer (we 
have verified also the case m = 2). I.e., T m is linearly 
uncorrelated with dT/dE and dT/dX. However, in gen- 
eral they are correlated; for instance we have cheked that 
(T(dT/dq) 2 ) ± 0. 
Following the same arguments we get 



dTdT\_ 1 
'dq'dq 7 / ~ Re 4 



E 



.6' = 1 



(('5'l)afc('5'r)a'6') 



diS^ab d(S 2 )a<b' 



(Si) al 



d(St) a , 



d(S* 2 )ab 

dq 



dq 



dq' 



(A4) 



where q and q' denote E or X independently. Also we 
have used that Si and S2 are unitary ((3 = 2) and sym- 
metric ((3 = 1) independent and equally distributed ran- 
dom matrices. When q = E and q' = X the result is 
((dT/dE)(dT/dX)) = because the partial derivative 
with respect to X is linear in the matrix Kj which has 
zero mean. Then, dT/dE and dT/dX are linearly un- 
correlated. Again, in general they are correlated. For 
instance, ((dT/dE) 2 (dT/dX) 2 ) ^0. 

Let us consider q = q' = X . Using the parametrization 
(|T5|l and averaging over Uj and Vj (j = 1, 2) for (3 = 2, 
or over Uj for (3=1, using the results of Ref & we arrive 
to 



((<5 , l)afc(<5'i)a'6' 



and 



(K'S b b +S b a S a b )/(N+l) (3 
S a a '8 b b '/N (3 



(3 = 1.2) 




= 1 
= 2 
(A5) 



(A6) 



for both (3 = 1, and 2, where we used (Qx aa ) = 
because Kj in the parametrization i|ltj|) has zero mean; 
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while 



for /3 = 1, 



/ d(S* 2 ) ab d(S 2 ) a , b , \ {5j5t + <%'<%') y /(C 2 s v 

N(N+1) 2-^\\^ x 'aa/ 



\ dX dX 
for (3=1, 



(A7) 



dX dX 



N 2 



J2((QxU) (A8) 



for = 2. Using Eqs. flU, JT^J) and JTSJl, Eqs. (£3 and 
can be written as 



/ d(S* 2 ) ab d(S 2 ) a , b . 
\ dX dX 



= 4 



iV(JV + l) 

AT iV 
E((^ 2 ) Q J+ E (QoaQ/3/3) 



(A9) 



for (3=1 and 

dX dX 



a'xb' N 



= 4- 



N 2 



E (Q««Q/3/3> (A10) 



for /3 = 2. 

We substitute Eqs. (|A5fl . (|A6j) . l|A9jl and (|A10I) into 
Eq. i|A4(l for q = q' = X . After we sum over a, a', b, b' 
the result is 



(dT/dX) 2 ) = 2 



N + l 



N N 
a=l a, P 

(All) 



1 N 

(dT/dxf) = - (Q aa Q/3 ) (Ai2) 



for 13 = 2. 

Now, the dialgonalization of Q, Eq. I|19(l . and the re- 
sults of Refi^i lead us to 



(dr/dxy 

for 0=1 and 



(A13) 



((9T/9A) 2 ) = (r?) + (N - 1) (nr 2 ) (A14) 
for /3 = 2. 

Finally, the average with respect to the t variables are 
done using Eq. I|14J) . When N = 1 only the first term of 
Eqs. (|A13|> and (|A14|I contribute. But, (rf) diverges for 
(3 = 1 and (3 = 2 and then ({dT/dX) 2 ). For AT > 1 we 
have 



2X _ / 2JV!/(iV-2)(JV + l)! /8 = 1 



and 



2AT(AT- 2)!/(AT+ 1)! ,9 = 2 



(A15) 



(riT 2 ) = (AT-l)!/(Ar + l)! = 1,2. (A16) 



The final result for {(dT/dX) 2 ) is given by Eqs. JH3 and 
for = 1 and /? = 2 respectively. 
In similar way we arrive to the result given by Eq. I|39() 
for the variance of dT/dE. 
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